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ABSTRACT 

Using two-dimensional velocity maps and /-band photometry, we have created mass models of 40 
spiral galaxies using the Milgrom relation (the basis of modified Newtonian dynamics, or MOND) 
to complement previous work. A Bayesian technique is employed to compare several different dark 
matter halo models to Milgrom and Newtonian models. Pseudo-isothermal dark matter halos provide 
the best statistical fits to the data in a majority of cases, while the Milgrom relation generally provides 
good fits as well. We also find that Milgrom models give mass-to-light ratios that roughly correlate 
with galaxy color, as predicted by stellar population models. A subsample of galaxies in the Hydra 
cluster follow a tight relation between mass-to-light and color, but one that is significantly different 
from relations found in previous studies. Ruling out the Milgrom relation with rotational kinematics is 
difficult due to systematic uncertainties in the observations as well as underlying model assumptions. 
We discuss in detail two galaxies for which the Milgrom relation appears to fail and find that relaxing 
the assumption of constant stellar mass-to-light ratio can maintain Milgrom models' viability. 

Subject headings: galaxies: fundamental parameters — galaxies: kinematics and dynamics — galaxies: 
dark matter — galaxies: stellar content 



1. INTRODUCTION 

The vast majority of spiral galaxy rotation curves 
are flatter than Newtonian pr edictions at la r ge dis- 
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1986f). Galaxies that do not show this behavior often 



show peculiar morphologies or perturbing interactions. 
This asymptotic flatness of rotation curves, which is not 
perf ect (for a more in-depth discussion of thi s point see, 
e.g.. iPersic et aZ.lll996l : ICatinella et aZ.I l2006h . disagrees 
with the Newtonian prediction in which the circular or- 
bital speed due to the gravity from a galaxy's luminous 
material should decrease at distances which enclose sub- 
stantial fractions of a galaxy's visible mass. This type 
of discrepancy is usually viewed as evidence for a non- 
luminous mass component, known as a dark matter halo. 
This picture of luminous matter embedded in dark mat- 
ter halos has become the paradigm of galact ic structure, 
and recent work involvin g galaxy clusters (|Clowe et all 
l2006t iBradac et aZll2006f ) has presented strong evidence 
that this accurately reflects astrophysical reality. Since a 
precise mass distribution of the visible matter is difficult 
to determine directly, the dark mat ter distribution is also 
not known. Previous stu dies (e.g.. iBlais-Ouellette et~al\ 
120011: iBarnes et q/.|[2004l . hereafter BSK) have found that 
the choice of dark matter distribution significantly im- 
pacts the derived mass of the luminous material. 
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A different phcnomcnological approach to rotation 
curve fitting is provided by the fitting form ula origi- 
nally proposed by Milgrom (|Milgrod[l983allbl lcll which 
he termed MOND, for "modified Newtonian dynamics." 
The original goal of MOND was to use only the lumi- 
nous matter to explain two widely observed properties 
of spiral galaxies; asymptotic flatness of galactic rota- 
tion curv es and the wel l -know n Tully-Fisher relation 
L ~ d /3 (jTullv fc Fisheri Il977f l. A tight relation be- 
tween a galaxy's luminosity L and the asymptotic value 
of the galaxy's circular speed v appears to be valid for 



well as barred and unbarred galaxies (ISpravberrv et al. 


1995 


;IZwaan et aZ.lll995 


:ICourteau & Rixlll999l;ICourteau 


2003 


). iMcGaugh et al. 


(2000f) has convincingly argued 



that luminosity is actually a stand-in for total bary- 
onic mass; i.e., the same Tully-Fisher relation holds 
for both low and high surface brightness spirals if bary- 
onic mass (stars and gas) is used in place of luminosity. 
When this bar yonic Tully-Fish er relation is examined, 
-^baryonic oc v A (|McGaughll2005f ) . These observations are 
built into the MOND formalism as follows. The total 
acceleration due to a gravitating object is given by the 
expression 

a/i(a/a ) = sl N , (1) 

which relates the total acceleration a to the Newtonian 
acceleration a^v through a function fi satisfying the lim- 
iting for ms ii — > 1 when a > oo and /z — > a/ao when 
(lMilgromlH983ah . The new universal constant 
ao = 1.2 x 10~ 10 km/s 2 is obtained from fits to mea- 
sured rotation curves and is the same for all galaxies. 
The exact form of the interpolating function is essentially 
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unconstrained by rotation curve data; throughout this 
paper we use n(x) = x/\fl + x 2 (|Begeman et aLlll991l) . 
It has long been realized that Equation |T]) violates con- 
servation of linear momentum; a more acceptable gravi- 
tational modification alte rs the Poisson equation for th e 
gravitational potential <p ([Bekenstein fe Milgrom|[l984f ). 



where p s is a characteristic density and r s is a scale 
length; pseudo-isothermal (PI) 
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This modification of the gravitational field explicitly con- 
serves energy and momentum, as it is based on a La- 
grangian formulation. However, in the limited case of 
spiral galaxies like those considered here, the force law 
modification of Equation (|T|) is a sufficiently accurate ap- 
proximation. 

In the inner regions of typical spiral galaxies, the char- 
acteristic acceleration scale a>ao and the gravitational 
force reduces to the usual Newtonian force. In the outer 
regions of spiral galaxies, a « ao (the so-called MOND 
regime) leading to the acceleration having a 1/r depen- 
dence, asymptotically flat rotation curves, and an en- 
closed mass varying like the flat rotation speed to the 
fourth power. Assuming a constant mass-to-light ratio 
T results in the Tully-Fisher relation. A number of pre- 
vious studies have utilized Equation (U) to describe ro- 
tation curves (lKentlll987l; iBegeman e^aZ-l 119911: ISanders I 
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modeling rotation curves with one free parameter (two if 
a bulge is present) rivals dark matter halo models that 
require at least two additional parameters. The empirical 
accomplishments suggest that regardless of the veracity 
of the hypothesis that Newtonian dynamics is modified, 
EquationQ] (the Milgrom relation) is a powerful and com- 
pact method for characterizing the rotational kinematics 
of spiral galaxies. 

In this work, we regard the Milgrom relation as an 
empirical fitting formula for rotation curve data. In par- 
ticular, we make no assumptions about properties of dark 
matter halos, or modifications to gravitation, from which 
th e Milgrom relatio n might arise. But as emphasized 
in iMcGaughl ([2004D , Equation flU is the effective force 
law for the visible matter in spiral galaxies, regardless of 
whether the force comes from a distribution of dark mat- 
ter or from MOND, which seems increasingl y unlikely in 
the fa c e of evidence li ke th at presented in Iciowe et al ' 
(120061) : IBradac et al\ (120061 ) (but see also lAngus et al 
2006). This paper remains agnostic about the ontolog- 
ical status of the Milgrom force law. However, we note 
that if galactic dark matter is taken as a certainty, then 
the success of the Milgrom relation in predicting spiral 
galaxy rotation curves reflects an underlying dynamical 
process in galaxy formation which is not well understood 
and is deserving of further scrutiny. In this work, we are 
interested in determining i) the extent to which this dy- 
namical relation holds and ii) specific situations in which 
the relation might not hold. 

This paper performs Milgrom fits to the same sample of 
spiral galaxies considered in BSK, and provides a statis- 
tical comparison with the dark matter hal o fits. The halo 
profile s considered in BSK are the familiar [Navarro et all 
(19961 hereafter NFW) 
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with central density po and core radius r c ; and a simple 
power law 

p(r)= P0 (j^j , (5) 

along with a Newtonian stars-only fit which is useful for 
comparing with both Milgrom and dark halo fits. 

Sample selection, photometric and kinematic data, and 
modeling procedures are discussed in § [2J Along with the 
results of Milgrom fits to this data 13. 1[) . we obtain the 
T-color relation for these models (§ 13. 2p . This sample of 
galaxies contains 8 Hydra cluster galaxies, and we focus 
on these galaxies to investigate the claim that Milgrom 
models produce T values that closely agree with popu- 
lation synthesis predictions and therefore provide likely 
estimates of luminous mass. We also compare the Mil- 
grom model fits to the dark halo fits from BSK, utilizing a 
Bayesian technique that quantitatively compares models 
with different numbers of free parameters. A descrip- 
tion of this technique along with model comparisons is 
given in Section 2] Section [S] discusses the impact of sys- 
tematic distance uncertainties. Finally, we discuss two 
anomalous galaxies in § [51 and consider the validity of 
these statistical goodness-of-fit measures in the presence 
of known systematic errors in our underlying assump- 
tions. Color gradients within individual galaxies and un- 
observed gas point to small variations in their mass-light 
ratios, at a level which naturally explains Milgrom fit 
discrepancies. 

2. DATA & MODELING 

BSK contains a detailed discussion of the data and 
fitting methods used in this analysis, and the interested 
reader is pointed there for a detailed description. Briefly, 
we use a sample of 40 high surface brightness galaxies 
with /-band images and Fabry-Perot two-dimensional ve- 
locity maps. These data have been provided by Povilas 
Palunas and Ted Williams and a thorough description of 
th e acquistion and reduction of the data can be found 
in iPalunas fc Williams! (|2000D . We decompose the I- 
band images into disk and bulge components when nec- 
essary. The disk is assumed to be radially symmetric 
and infmitcsimally thin and the bulge is assumed to be 
spheroidal. 

A mass model is determined once the mass-light ratios 
for the disk (Tdj) and the bulge (Tg/) are s pecified. 
Desp ite evidence for gradients in T factors (e.g., Ide Jond 
we first consider T values to be spatially constant; 
we relax this assumption in § 16.31 The / subscript indi- 
cates mass-light ratios determined from /-band photom- 
etry. We neglect any contribution from gas, as we have 
no information about the radial distribution of gas in 
the galaxies in our sample; however, the gas in late-type 
galaxies like those in our sample is likely to contribute 
< 10% of the ci rcular speed at any radius (|Broeilslll992l : 
IVerheiienI Il997t ) . Neglecting it will cause a slight sys- 
tematic overestimate of Td,i and can lead to models 
mildly underpredicting circular speeds in the outer parts 
of galaxies. It could be argued that omitting the gas 
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contribution will cause our stars-only Milgrom models 
to describe the data worse than other models simply be- 
cause we are not including all relevant information; i.e., 
our Milgrom models are born at a disadvantage. How- 
ever, our Milgrom models do not show any systematic 
trends indicating that large amounts of gas would im- 
prove the models; the Milgrom models do not systemati- 
cally underpredict the rotation speeds at the edge of the 
kinematic data. We discuss the impact of this omission 
in more detail for two interesting galaxies in § [SJ 

The disk-bulge mass model determines the Newtonian 
circular speed v c (6) as a function of angular separation 
from the center of the galaxy. This is converted to a cir- 
cular speed v c (r) as a function of distance from the center 
of the galaxy, using the angular diameter distance to the 
galaxy given by the Hubble distance (systematic errors 
associated with this simple distance assumption are dis- 
cussed in Sec. 15. ip . Equation [T] converts this Newtonian 
circular velocity to a Milgrom-relation circular velocity. 
The Newtonian disk-bulge mass model is also combined 
with the standard dark matter halo models (Eqs. OQJ 
and [5]) to give the total rotation curves in these cases. 

The rotation curve inferred from the model mass dis- 
tribution is compared with the actual galaxy rotation as 
measured by the Fabry-Perot maps. Throughout this 
paper, we fix the position angle <f>Q, inclination angle i, 
and systemic velocity v s of each galaxy at the values 
deduced from earlier fit s of axisymmetric flow pa tterns 
to our kinematic data ([Barnes fc Sellwoodll2003h . We 
have found the values of these parameters that mini- 
mize x 2 between the two-dimensional velocity map and 
a planar axisymmetric flow pattern. The fit therefore 
yields a set of mean orbital speeds, and their uncertain- 
ties, at equally-spaced radii which we show as the "data" 
in subsequent figures illustrating rotation curves. Since 
we fit axisymmetric mass models, we could, in principle, 
fit only this onc-dimentional rotation curve, which rep- 
resents the speed averaged in annular bins, provided the 
weights are computed correctly. However, the full two- 
dimensional maps have several important advantages, in- 
cluding clear identification of outlying pixels and identi- 
fication of non-axisymmetric velocity patterns (usually 
bars or spirals) which cause systematic errors in rotation 
curves obtained from one-dimensional velocity measure- 
ments. Of course, identifying non-axisymmetric patterns 
is only half the battle. We are looking towards collecting 
data for a sample chosen with minimal spiral and/or bar 
structures that will involve fewer systematic effects and 
will allow better model discrimination. 

3. MILGROM FITS 

The mass model incorporating disk and bulge compo- 
nents gives a Newtonian circular velocity v c (r) and accel- 
eration ajsr(r). Equation Q] then gives a total acceleration 
a(r); equating these accelerations with centripetal accel- 
erations give us the Milgrom circular speeds. Note that 
for the galaxies considered here, the simple acceleration 
approximation in Eq. [T] gives essentially the same an- 
swer as the physically more realistic Poisson-like Eq. [2] 
which is numerically more challenging to solve. Using 
previou sly determined projection angles and systemic ve- 
locities ([Barnes fc Sellwoodll2003h . which are held fixed, 
we project a one-dimensional speed profile into a two- 
dimensional velocity map. This model velocity map is 



compared to the data velocity map, pixel by pixel. 

A small number of pixels exhibit what are clearly 
anomalous inferred velocities; these are almost all iso- 
lated pixels in outer regions of galaxies, where the ve- 
locity is much different than the smooth velocity field 
defined by the surrounding pixels. These points arc due 
to imperfect filtering of cosmic rays and unmodcllcd line 
emission in the data analysis pipeline, which spuriously 
fits the pixels as having high velocity and low dispersion. 
To reduce the influence of these non-normal errors, we 
find the best-fit values for Td,i and Tgj by min imizing 
the Tukey biweight parameter ([Press 
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where z, ; = (Dj — Mj)/<7j, Di and Mj are the respective 
measured and mass-model velocities in pixel i, Gi is the 
measurement standard error on D„ c is a constant, and 
N is the number of data points. Note that in the limit 
c — > 00, this reduces to the usual expression for reduced 

chi-squared xl = U N T,Zii D i ~ A*i) 2 /°i- While data 
values that differ from the model prediction by < cai/2 
contribute to xl almost exactly as for the conventional 
Xr, the contribution from data values that are farther 
from the model prediction by co~i does not change as the 
parameters are adjusted, and outlying data therefore do 
not in fluence the resulting parameter values. iPress et all 
(|1992f ) recommend c = 6, but non-circular flow patterns 
in our galaxies suggest that a larger value would be more 
appropriate; after some experimentation we find c = 10 
eliminates the influence of the extreme outlying veloc- 
ities while allowing the great majority of the pixels to 
contribute to the fit with almost full weight. Only a 
handful of pixels per galaxy image arc ignored using this 
procedure. 

3.1. Parameter Values and Rotation Curves 

The best-fit \b> ^d,i and Tb,i parameters, statisti- 
cal uncertainties (a), and systematic uncertainties (s) 
for the Milgrom mass models are given in Tabic [1] We 
also include the ratio of the distance to the edge of the 
kinematic data i?k,max to the disk scalelength rj for 
each galaxy. On average, the kinematic data extends 
to f=a 4.5rd and covers the visible disk. The uncertainty 
calculations are discussed in the next paragraph. The av- 
erage Td,i is 1.9 and the average Tb,i is 2.9, with wide 
variation in both values. Generally, the Milgrom model 
values are lower than the Newtonian stars-only values; 
the ratio of T between the Milgrom models described 
here and the Newtonian stars-only models of BSK are 
0.6 and 0.5 for the disk and bulge components, respec- 
tively. 

Statistical uncertainties have been calculated using 
Markov Chain Monte Carlo (MCMC) techniq ues (e.g., 
iChristensen et al\ |2001| ; iKosowskv et al\ l2002f ) to map 
out allowed regions of parameter space. Projections of 
the boundary which encloses the most likely 68% of the 
models then provide the statistical uncertainties. The 
average statistical uncertainties are 0.08 for T^j and 
0.3 for Tb.i- The quoted systematic uncertainties in 
the fitted Tj values are due to the uncertainties in the 
galaxy's inclination and position angle. A previous pa- 
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per (|Barnes fc Sellwoodll2003f ) estimated systematic un- 
certainties in inclination and position angles caused by 
spirals and other nonaxisymmetric structure in the galax- 
ies. The best-fit value plus the estimated high and low 
limits for each of the two projection angles yield nine 
different combinations of inclination and position angle 
which give nine separate estimates for Tj values. We 
take the systematic uncertainty in a T/ value to be half 
the range of these nine Tj values. As in BSK, the sys- 
tematic uncertainties dwarf the statistical ones with an 
average sr D 1 of 0.5 and an average sr B , of 0.8. 

Figure [T] illustrates the rotation curves of six galax- 
ies (assuming fixed Hubble distances). The best-fit Mil- 
grom models (solid lines) are shown along with the data 
points derived from velocity maps. For comparison, wc 
also show the PI halo models (dot-dashed lines) from 
BSK. The top two panels are for galaxies in which the 
Milgrom fit is favored over the various BSK halo models. 
The middle set of panels show fits to galaxies with non- 
axisymmetric structure in which the Milgrom fit is less 
favored than the other BSK models. The bottom pan- 
els show curves for axisymmctric galaxies in which the 
Milgrom fit is less favored than the BSK models. These 
galaxies will be discussed in more detail in §[51 One of the 
most striking things about this figure is the close similar- 
ity bewteen the Milgrom and halo curves; the differences 
in circular speeds between the models shown are rarely 
much larger than 10 km/s. Large differences in likeli- 
hoods between the models are due to small variations in 
model velocities over very large numbers of pixels in the 
two-dimensional velocity maps. 

3.2. Color-T relations 

Stellar population synthesis models predict correla- 
tions between galaxy color and Ypj value, with red 
galax ies having larger Tp / values than blue galaxies 
fe.o.. IBell fc de Jong|[200llBell et aLll200l . As in BSK, 
we have created galaxy colors obtained from the ESO- 
LV aperture magnitudes, as recorded in the NASA/IPAC 
Extragalactic Database (NED). In principle, we should 
compare the Tjj j with the color of the disk alone, but 
NED gives magnitudes of the disk and bulge combined. 
The bulge light in our mostly late- type galaxies is a small 
fraction, typically < 10%, of the total in the J-band, so 
any correction for the bulge to the total color is likely 
to be small. We have corrected the NE D magnitudes 
for ga lactic extinction using values from Schlcgcl et al. 
(11998T) and for i ntern al extinction using the prescriptions 
of iTullv et all (|1998l ). We estimate errors in our colors 
to be around 0.15 mag. 

In BSK, we found that the model T values for the 
galaxies in this sample did not have any trend, and in- 
stead follow a flat distribution, albeit with large scatter. 
We present here the color- T^ j relation for the Milgrom 
models for these same galaxies. Figure [2] shows the T o j 
value of each galaxy plotted versus its dereddened color 
for Newtonian stars-only models (top frame) and Mil- 
grom models (bottom frame) . The error bars reflect the 
statistical and systematic uncertainties added in quadra- 
ture. Galaxies marked by squares have no bulge compo- 
nents. The s olid lines in each f rame represent the correla- 
tion given in iBell et all (|2003h . The dashed lines are the 
best-fit linear fits to the data. Contrary to the findings 
of BSK, we see a trend in the Milgrom color-T plane sim- 



ilar to that expected by stellar population models. The 
slopes of these lines show some discrepancy but there is 
still a large scatter in our points. 

This sample of 40 galaxies contains 8 which belong to 
the Hydra cluster. The color-T p.i plot of these gal axies 
is shown in Figure[3] As before, the lBell et all ()2003l ) cor- 
relation is shown as a solid line, while the best linear fit to 
the data (marked by diamonds) is the dashed line. This 
linear fit appears to accurately describe the correlation 
present in the data with little scatter. The large offset 
between the two lines is surprising since lMcGaughl (I2004D 
found that MOND values matched the IBell et all 
(|2003| ) line remarkably well. There are a number of pos- 
sible explanatio ns of the d i screpa ncy between the Hydra 
galaxies and the lBell et all (|2003h correlation. One possi- 
bility is that the offset is due to errors in the adopted dis- 
tances to these galaxies (see § [5] for a detailed discussion 
of th e impact of di s tance on MOND T values). How- 
ever, iHudson et all (|l997f l give a Fundamental Plane- 
based, Malmquist bias corrected distance to the Hydra 
cluster of 52.6 ± 2.5 Mpc and the galaxies in this plot 
have adopted distances that average to 52.6 ± 4.2 Mpc. 
Metallicities and initial mass functions of galaxies have 
also been shown to strongly impa ct the slope and inter - 
cept of the color-Tuj correlation (|Bell fc de Jondl200lh ; 
these possibilities would require further observations to 
explore in more detail. 

4. MODEL COMPARISONS 

Previous studies have rarely attempted direct compar- 
isons between Milgrom fits and halo model fits. Here 
we compare the Bayesian probabilities of halo models to 
Milgrom models, taking into account the differing num- 
ber of free parameters in the models. A brief summary 
of the necessary Bayesian probability theory is given in 
the Appendix. The model comparisons here are based 
on statistical fits only. Evaluation of some systematic 
sources of error and their potential impact on straight- 
forward statistical model comparison arc considered in 
the following sections. 

We are considering several different models with dif- 
ferent numbers of parameters and we would like to know 
which provides the most economical description of the 
data. Specifically, we have compared our Milgrom fits to 
various halo fits (PI, NFW, and power law halos), along 
with Newtonian stars-only fits for illustration (for a full 
discussion of the various halo models, see BSK). We have 
used the following parameter ranges for (uniform) prior 
probability values: (PI) Ap = 2 M /pc 3 , Ar c = 20 
kpc; (NFW) Ac = 20, Au 20 o = 300 km/s; (power law) 
A^o = 300 km/s, Aa = 2. Altering the priors will change 
a given model's probability, but other sensible choices do 
not change our basic results. Tabic [5] lists the logarithms 
of the probability ratios A of the BSK models to Milgrom 
fits. A values < indicate that the Milgrom model is pre- 
ferred. The large ratios between some halo and Milgrom 
models are due to the large numbers of data points in 
the two-dimensional velocity data, resulting in very small 
statistical error bars and large likelihood ratios between 
models with similar rotation curves. A simpler model 
compar ison test bas ed on Bayesian information criteria 
(BIC) (|Liddlell200l agrees with the results we present 
here. 

Surprisingly, the Newtonian stars-only model is 
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strongly favored over Milgrom for nine galaxies, and the 
models are comparable for two others (ESO 317g41 and 
445g39). In fact, the Newtonian stars-only model is pre- 
ferred over all other models for ESO 438g08 and 382g06. 
Of these 11 galaxies, 5 have substantial bars, 3 have no- 
ticeable spirals, and 1 has a close, possibly interacting, 
companion (ESO 438g08). Since it is possible that these 
circumstances introduce a bias into the model compari- 
son, we do not focus our attention on these galaxies. The 
two remaining galaxies (ESO 323g42 and 382g06) appear 
axisymmetric (no large bars or strong spiral structure). 
These galaxies are discussed in detail in § [6] 

For three galaxies, Milgrom models are preferred or 
are comparable to all other models: ESO 322g76, ESO 
445gl5, and ESO 501g68. ESO 445gl5 and ESO 501g68 
have bars while ESO 322g76 has moderately strong spi- 
ral structure. All 3 have velocity maps that extend at 
least 2.5 disk scalelengths and their rotation curves have 
flattened by the last data point. Of the 40 galaxies an- 
alyzed, these are the only galaxies for which Milgrom is 
preferred over the PI halos. Milgrom models are pre- 
ferred over NFW halos in 7 cases. These results suggest 
that, in general, the simplicity (fewer free parameters) 
of Milgrom models does not statistically outweigh the 
ability of halo models to provide better fits. 

Due to its poor performance relative to the halo mod- 
els, it is tempting to conclude immediately that the Mil- 
grom relation describes galaxy rotation curves only by 
accident, undermining the strongest reason for consider- 
ing MOND as a modification of gravity. However, this 
conclusion would be premature without detailed consid- 
eration of systematic uncertainties. It is dfhcult to find 
cases where the poorer fit of the Milgrom model cannot 
be accounted for by reasonable known systematic errors. 

Milgrom models are highly constrained, and we have 
made several implicit assumptions which have substan- 
tial effects on the rotation curves: 1) all stellar motions 
are circular; 2) the visible stellar light accounts for all 
of the mass in the galactic disk (neglecting, in particu- 
lar, any gas, which is known to be non-negligible in some 
galaxies); 3) the mass-light ratio is constant throughout 
the galaxy; 4) the adopted inclination angle is correct; 
and 5) the galaxy distance is given by the Hubble dis- 
tance. Not one of these assumptions are exactly true, 
and in many galaxies these assumptions are known to 
be incorrect at levels which would change the rotation 
curve by an amount significantly larger than the statis- 
tical errors. The sample considered here, in particular, 
was chosen to probe large-scale velocities in the local uni- 
verse, not because of its particular suitability for galaxy 
dynamics; we have analyzed it because of its large sample 
of high-quality two-dimensional velocity measurements. 
Milgrom models have essentially no ability to compen- 
sate for any failure of these assumptions as they are de- 
termined completely by the distribution of visible matter. 
Halo models, by contrast, have two degrees of freedom 
which are unconnected to the distribution of visible mat- 
ter and otherwise unconstrained by direct observation, 
and can therefore compensate to some degree for any in- 
correct assumptions about the visible matter. In other 
words, beyond the fact that halo models will generically 
provide better fits because they have more parameters, 
they are also more flexible in compensating for system- 
atic errors, which is not reflected in purely statistical 



comparisons. 

In the remainder of this paper, we consider in more de- 
tail some of the important systematic errors which might 
affect the analysis presented here. We then focus on two 
sample cases for which Newtonian fits are statistically 
superior to either Milgrom or dark halo fits, and discuss 
which possible explanations for the discrepancy are rea- 
sonable. 

5. SYSTEMATIC EFFECTS 
5.1. Distances 

A substantial source of uncertainty is the adopted dis- 
tance to a galaxy. T values arc inversely proportional 
to distance in both Newtonian and Milgrom mass mod- 
els. Unfortunately, a significant fraction (« 50%) of the 
galaxies in this sample were originally observed to study 
peculiar velocities in the region toward the "Great At- 
tractor" and may therefore have larger-th an-usual pecu- 
liar ve locities and distance uncertainties. iBothun et all 
(|1992t ) used the Tully-Fisher relation to estimate dis- 
tances to this subsample of galaxies, finding differences 
between redshift distance and Tully-Fisher distance of 
typically about 25%. While for Newtonian (stars-only 
and halo) fits this simply introduces a scale uncertainty 
in T values, Milgrom fits are subject to a further dis- 
tance dependence, since they depend on the dimensional 
parameter ag. Changing the distance to a galaxy changes 
the acceleration scale at a given angular separation from 
the galaxy center, and thus changes the inferred radius 
at which the non-Newtonian behavior sets in. For exam- 
ple, suppose the distance to a galaxy has been adopted at 
half its actual distance. Given angular separations now 
correspond to physical distances that are half the ac- 
tual values. Since the accelerations are inversely related 
to physical distances, the accelerations increase. These 
increases then push the location where acceleration be- 
comes comparable to ao farther out in the disk. Mov- 
ing a galaxy closer forces the Milgrom fit to be closer 
to Newtonian, while adopting a distance larger than the 
actual value increases the discrepancy between Milgrom 
and Newtonian fits. 

While Milgrom fits are fairly sensitive to the 
adopted distance scale, this fact can be exploited as 
a method for determining distances. As related in 
ISanders fc McGaughl ([2002) , changing the adopted dis- 
tance to a galaxy until a satisfactory Milgrom fit is ob- 
tained gives an independent distance measurement that 
is often in accordance with Tully-Fisher distances (which 
is not surprising, as the Tully-Fisher relation is implied 
by the Milgrom relation) . We have investigated this pos- 
sibility by including the distance to galaxies as another 
free parameter in Milgrom fits. For our data set, most fits 
of this kind are unstable, leading to models with galax- 
ies unreasonably nearby or far away. Part of the reason 
for this behavior is due to bars and spirals. Nonaxisym- 
metric features tend to lower rotational velocities, mak- 
ing axisymmetric models overshoot the data, especially 
in central regions. In such a case, the fitting routine 
can better match the data by moving the galaxy far- 
ther away, thereby lowering the T values, making the 
baryonic component insignificant, and allowing the non- 
Newtonian MOND behavior to dominate. This is sim- 
ilar to the behavior of some galaxies that are best fit 
with a halo alone, as noted in BSK. One of our galax- 
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ies (ESO 438g08) has a declining rotation curve that is 
best fit by the Newtonian stars-only model and returns 
a terrible fixed-distance MOND fit. Once the distance 
is allowed to change, the galaxy is moved unbclicvcably 
closer so that the fit can become more Newtonian. If this 
galaxy's kinematics can be shown to be due only to the 
matter within the galaxy, it would pose quite a problem 
for MOND, as well as for any commonly considered dark 
matter halo. However, the galaxy has a nearby, possibly 
interacting, companion; while we cannot definitively say 
that tidal interaction is causing the declining rotation 
curve, neither can we rule it out. 

In an attempt to remove the effects of aggravating cen- 
tral asymmetries, we have excised pixels from the centers 
of the galaxy images and re-fit the remaining data with 
a Milgrom model that has the distance as an additional 
free parameter. If the free parameter distance varies by 
more than a factor of 3 from the Hubble distance we 
used in the fixed distance fitting, we consider the dis- 
tance fit unstable. We fit each galaxy 10 times, each 
time increasing the number of excised pixels. The first 
attempt utilizes all the pixels and each successive trial 
removes an additional 5% of the total radial data range 
from the central regions. Fifteen galaxies cannot find ac- 
ceptable distance fits for any amount of excised pixels. 
Fourteen galaxies have acceptable fits using all the data. 
Twenty-two galaxies have at least three consecutive fits 
with acceptable distances; however, the distances that 
these fits return never agree with each other to better 
than 30% of the Hubble distance. So while the adopted 
Hubble distance necessarily introduces a systematic un- 
certainty that predominantly affects the Milgrom fits, we 
have not been able to devise a stable method of deter- 
mining more accurate values directly from the dynamics 
of the galaxies. Rotation curve data extending further 
into the flat portion of the rotation curve force a stable 
distance determination from Milgrom fits; this is clearly 
an important limitation of the data set considered here. 

5.2. Projection Angles 

Systematic uncertainties in projection angles substan- 
tially affect derived T values (see the discussion in § 13. ip . 
However, the estimates of systematic uncertainties in the 
projection angles appear to be reasonably behaved. Our 
sample contains 14 galaxies for which the uncertainty in 
T from the estimated systematic uncertainty in projec- 
tion angles is less than twice the statistical uncertainty 
(see Table 1); in most of these cases the systematic error 
is less than the statistical error. A cursory examination 
of the model fits for these 14 galaxies shows that the var- 
ious halo models or the Milgrom model are preferred as 
the best-fitting model in roughly the same fractions as 
for the full set of galaxies. More importantly, the galax- 
ies of this subsample tend to have |A| < 10, while the 
full sample contains several cases with |A| > 100. These 
large values are likely a signal that the straightforward 
statistical comparison may be suspect, and that system- 
atic errors arising from projection angles or other sources 
are contributing to the strong rejection of the Milgrom 
model. 

5.3. Nonaxisymmetry 

Many of the galaxies in the sample have bar and/or 
spiral features or lopsided shapes. In these cases, the 



assumption of axisymmetry is manifestly violated, al- 
though if the irregularities are weak (as is generally 
the case) the effect on the rotation curve will be small. 
Departures from regular elliptical (projected) isophotal 
shapes can indicate recent merging or tidal interactions 
which substantially distort the velocity field. These dis- 
tortions are evident in maps of the differences between 
circ ular orbit models and kine matic data {e.g., see Figure 
5 in iBarnes fc S cllwood 2003]). The effect of departures 
from axisymmetry requires painstaking modelling (e.g., 
IWeiner et all |2001D . but incisive comparisons of halo 
models with Milgrom fits should be limited to galaxies 
with no gross departures from axisymmetry. Our sample 
contains roughly equal numbers of galaxies with strong 
(18) and weak (22) spiral structure. Seven galaxies have 
some evidence for a bar; 5 are strong bars in galaxies 
with strong spiral structure, and 2 are short/weak bars in 
galaxies with weak spiral structure. The specific galaxies 
discussed in the following section have only weak spiral 
structure. 

6. ESO 323G42 & ESO 382G06 

MOND is vastly more falsifiable than dark matter hy- 
potheses. Since the fitting parameters are few and un- 
derstood more easily than the usual halo parameters, any 
one galaxy that can not be fit by the Milgrom relation 
could, in principle, refute the idea. Most galaxies in our 
sample have better halo fits than Milgrom fits; however, 
as discussed in § 15. li these relative likelihoods may be 
the result of systematic errors. 

A better way to rule out MOND is simply to find galax- 
ies with poor Milgrom fits in an absolute sense. This 
criterion also is hard to define as measures of fit qual- 
ity are also affected by systematic errors. One simple 
way to judge the quality of a Milgrom fit is in a rela- 
tive sense, by comparing with the Newtonian stars-only 
model. If the Milgrom relation provides an accurate de- 
scription of galaxies, then Newtonian fits should never 
be better than Milgrom fits (they will be comparable if 
all of the measured velocities are close enough to the 
galaxy center that they are still in the Newtonian grav- 
itational limit of the Milgrom force law). Many of the 
systematic effects that the halo fits can compensate for 
will affect both Milgrom fits and Newtonian fits equally, 
since both have little model freedom to compensate for 
systematic errors. It is difficult to see, for example, how 
an unaccounted-for distribution of gas could improve the 
fit of a Newtonian model relative to a Milgrom model; 
it should affect both equally and in the same direction 
relative to the measured velocities. 

Finding galaxies whose Newtonian fits are clearly su- 
perior to Milgrom fits can kill MOND as a fundamen- 
tal theory. (Note that galaxies whose Newtonian fits 
are superior to dark halo fits also pose severe problems 
for dark matter models.) Here we focus on two galax- 
ies, ESO 323g42 and ESO 382g06, whose Milgrom fits 
are worse than their Newtonian stars-only fits. Their 
rotation curves and several model curves are shown in 
Fig-HJa) and (b), respectively. 

In practice, the path to ruling out MOND is not so 
clear as systematic uncertainties complicate the situa- 
tion. For these two galaxies, the residuals from our kine- 
matic fitting do not show organized patterns associated 
with strong bars and/or spirals, and we believe these 
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galaxies are quite axisymmetric. Any spiral structure 
present likely perturbs circular speeds by < 10 km/s, 
and the resulting small rotation curve features would 
be similar in both Newtonian and Milgrom fits. With- 
out detailed information about the gas distribution, it 
is certainly possible that these Milgrom models are per- 
forming poorly simply because they have incorrect in- 
put. Further, the observed discrepancies between the 
Milgrom models and the data have magnitudes ( < 10%) 
that c ould possibly be lifted by including the gas (|Broeilsl 
Il992t IVerheiien|[l997t ). We note that this is a question 
that can be addressed with further observations. 

6.1. ESO 323g42 

Figure BJa) shows that the Milgrom relation fit (solid 
line) is clearly a poorer approximation to the data than 
either the Newtonian stars-only fit (dotted line) or the 
PI fit (dot-dashed line and the best halo fit) over the 
entire radial range of the data, w 3.4 disk scalelengths 
Rd (the NFW curve has been omitted for clarity but is 
very similar to the PI curve). The systematic projection 
angle uncertainties for this disk are quite small (« 2°) 
and likely do not impact the fit strongly. Unfortunately, 
this galaxy is also one of the cases where no stable fit 
could be found when the galaxy distance is left as a free 
parameter. The Hubble distance may not be the correct 
value, but there is no reliable evidence to change it. The 
discrepancy between the Milgrom curve and the data is 
most evident where the model undcrprcdicts v c near 5 
kpc and overpredicts v c beyond 12 kpc. A reasonable 
gas distribution could rectify these discrepancies. How- 
ever, one problem with adding significant amounts of gas 
to the model is that the model T o.i value necessarily has 
to decrease. The Ydj value in Table Q] is 0.78, already 
somewhat lower than the average value. Looking at the 
PI model curve, it is apparent that departures from ax- 
isymmetry arc quite small, < 5 km/s. Any mild spiral 
structure that may be present is not greatly influencing 
the model. 

6.2. ESO 382g06 

Unlike the case of ESO 323g42, the Milgrom fit for 
ESO 382g06 basically agrees with the both PI and New- 
tonian stars-only fits everywhere (out to the last point 
at w 4.3i?d), except in the very central region of the 
galaxy. In that region, no model explains the data well. 
One possible explanation is that, despite lacking a clear 
signature in either the photometric or kinematic resid- 
uals, the center of ESO 382g06 harbors some kind of 
non-axisymmetric structure. The rotation curve of ESO 
382g06 decreases slightly towards the edge of the opti- 
cal d isk, which is not so unco mmon for late type spi- 
rals (jBeauvais fc Bothunl I2001D . Such behavior in the 
low-acceleration regime causes the Milgrom relation fit 
to underpcrform the halo models, since it unavoidably 
gives asymptotically flat rotation curves. Note that this 
discrepancy is directly tied to the shape of the rotation 
curve, and will persist even including an undetected gas 
component. The projection angle uncertainties for this 
ga laxy are about average for this sample (as determined 
in iBarnes fc Sellwoodl [2003) . ~ 5°. However, it is un- 
likely that changing the projection angles will change the 
slope of the outer part of the rotation curve. Following 



the procedure described in § [5j differing amounts of cen- 
tral pixels have been removed from the fit, and in three 
trials stable distance values have been found. Unfortu- 
nately, these values do not agree with each other very 
well and the fits give unphysically small amounts of mass 
to the disk. Consequently, we again find that changing 
the adopted distance to the galaxy will not improve the 
Milgrom fit. 

6.3. Variable Ydj Fits 

In Section El we have discussed some of the observa- 
tional systematic uncertainties that could adversely in- 
fluence the Milgrom fits. Here we consider the impact 
of a systematic modeling error due to the fundamental 
assumption of a constant mass- light ratio T. Milgrom 
fits will be particularly susceptible to this type of error 
because of the deterministic dependence of the rotation 
curve on the baryonic mass. Halo models of any kind 
generally have sufficient freedom to contort themselves 
so that mild T gradients leave the quality of a halo fit 
unchanged. 

We have performed Milgrom fits for three galaxies 
(ESO 322g76, ESO 323g42, and ESO 382g06) allowing 
the Ydj values to vary radially; Y bj values continue to 
be treated as spatially constant. Specifically, we assume 
a linear relation, 

Tflj(r) = T_d, + mr, (7) 

where T^o is the central disk T value and to is the slope 
of the relation in (M0/L Q )/kpc. The rotation curves for 
the variable Ydj mass models are compared to the data 
points, the constant T curves, and the PI halo curves 
in Figure El The solid lines correspond to the variable 
Yd j fits, the dot-dashed lines are the PI halo fits, and 
the dotted lines are the constant T fits. 

The fit for ESO 322g76 (Figure Ek) serves as a test 
case since the constant T Milgrom fit is comparable or 
preferred to all the other models. We find T^o = 1.30, 
Y b,i = 1.75, and to = 3.5 x 10~ 3 . Since the slope is small 
and the T values are nearly unchanged from their original 
values (Yd, i — 1.33, Yb,i = 1-73), it is not surprising 
that the fits arc equally good (for both fits, \l = 2.05). 
For this galaxy, with an acceptable Milgrom model with 
constant T, the additional freedom does little to change 
the previous fit. 

This is not the case for the two galaxies that appear to 
have problematic Milgrom fits. Allowing for Ydj gra- 
dients, ESO 323g42 (Figure Eb) is best fit with the pa- 
rameters T Dj0 = Y B ,i = 1.16, to = -4.0 x 10~ 2 . The 
resulting xl value is 1.54, roughly 90% of the constant 
Ydj value, making it comparable to the constant T halo 
and Newtonian fits. The change in Ydj value from the 
center to the edge of the kinematic data is 0.29, much 
larger than cither uncertainty given in Table [1] This 
suggests that including the T gradient is truly improv- 
ing the fit for this galaxy. The best-fit parameters for 
ESO 382g06 (Figure Efc) arc Y D ,o = Y B j = 3.81 and 
to = —2.5 x 10 _1 . This case has a slight decrease in \\ 
value, from 2.48 with constant Ydj values to 2.46. In 
this case, the Ydj value changes from center to edge by 
1.2, smaller than just the statistical uncertainty alone in 
Table [TJ While the quality of the fit has definitely im- 
proved, it is not as clear in this case that allowing a T 
gradient has significantly changed the model. 
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Qualitatively, the Td,i slopes for these two galaxies 
act in the direction expected from observations of color 
gradients, with the oute rmost regions having lower T d.i 
values than the centers (|de Jondll996f ). Multi-band pho- 
tometry would allow determination of the actual color 
gradient, and whether this gradient is consistent with 
the best fit value. 

ESO 323g42 and ESO 382g06 exemplify the difficulties 
associated with falsifying the Milgrom relation. While 
there are certainly hints of problems with the zeroth- 
order models, plausible systematic effects (distance, gas 
mass, T gradients) must be exhausted before a definitive 
judgement can be handed down. For Milgrom fits with 
non-constant mass-light ratios, further observations can 
provide consistency checks, unlike halo models where all 
of our ignorance is folded into halo parameters and inde- 
pendent direct measurements of individual halo shapes 
are generally not possible. 

7. SUMMARY & CONCLUSIONS 

We have investigated MOND mass models of 40 high 
surface surface brightness spiral galaxies. Using two- 
dimensional Fabry-Pcrot velocity maps and /-band pho- 
tometric images, best-fit bulge and disk T values have 
been determined using the Milgrom relation (Equa- 
tion [T]). As in our previous work using dark matter halo 
models (|Barnes "eTZI[200l BSK), we find that the sys- 
tematic uncertainties due to projection angle uncertain- 
ties are typically much larger than the statistical uncer- 
tainties, ft should be kept in mind that large system- 
atic uncertainties could mask differences between vari- 
ous models and make any statistical comparison ques- 
tionable. 

We find that Milgrom model T o j values correlate with 
B — R galaxy col or (Figure Wp), in line with populatio n 
synthesis models (jBell fc de Jondl2001l: iBell et q/.ll2003D . 
While the correlation does not exactly match the syn- 
thesis model trend and contains much scatter, it sharply 
contrasts with the results of BSK where no trend was 
found (Figure [2jt). For a subsample of 8 galaxies that 
lie in the Hydra cluster, we find that they define a fairly 
tight Tuj-color correlation (Figure [3]). Again, the cor- 
relation does not agree with the trends from population 
synthesis models in either slope or intercep t. This con- 
trasts with the findings in McGaugh (2004) which show 
excellent agreement between Milgrom model and pop- 
ulation synthesis T-color r e lation s. We note that the 
galaxies used in iMcGaughl (|2004h are the same galax- 
ies ijosed_to_calil2rate_the population synthesis prediction 
in IBell fc de Jond (|2001h . We feel confident that the dis- 
agreement we see can not stem from distance errors, since 
the mean of the distances to these galaxies is the same 
as an inde pendent determinat ion using the Fundamen- 
tal Plane (|Hudson et aZ.I fl997t ) . It could be argued that 
bulge contamination skews our results, but these are all 
late-type spirals with very little bulge. It is difficult to 
imagine that these bulges could cause such a radical shift 
in the correlation. Recognizing that the common view is 
that IMFs are universal, one pos sible explanation made 
plain by the model calculations of lBell fc de Jonel (|200I[ ) 
is that there are different IMFs for these two sets of galax- 
ies. With only 8 galaxies, the evidence is weak to chal- 
lenge the standard view. 

Utilizing a Bayesian model comparison technique (§[4] 



and § [SJ , we have compared Milgrom models to Newto- 
nian stars-only, pseudo-isothermal halo, Navarro-Frenk- 
White halo, and power law halo models. The results are 
presented in Table [2] Overall, halo models arc statisti- 
cally favored over Milgrom models. At the same time, 
Milgrom models are generally preferred to Newtonian 
stars-only models. Unfortunately, this kind of analysis 
takes statistical uncertainties only into account. The 
quality of the model fits is also influenced by systematic 
observational uncertainties (distance, inclination, etc.; 
see § [5]). The Milgrom models are also subject to sys- 
tematic modeling effects (like assuming constant mass- 
to-light ratios) that are more easily absorbed in halo 
models. 

We have also discussed two specific galaxies which ap- 
pear to present a problem for the Milgrom relation (§|HJ). 
Our investigation suggests that nonaxisymmetry, projec- 
tion angle uncertainties, or distance errors are not re- 
sponsible for the Newtonian stars-only fits being better 
than the Milgrom fits. These galaxies provide excel- 
lent tests of the Milgrom relation since in this picture 
no galaxy should show Newtonian behavior in the low- 
acceleration (large radius) regime. 

With an eye towards explaining away the seemingly 
Newtonian behavior in these two galaxies, we have per- 
formed additional fits allowing the T of the disk to vary 
linearly with radius. We find that for the galaxy ESO 
322g76, which is statistically well explained by the Mil- 
grom relation, the best-fit parameters of a constant T^^ 
are virtually unchanged by allowing Td,i to vary. How- 
ever, for the two problem galaxies ESO 323g42 and ESO 
382g06 we find that 1) the best-fit parameters of the 
varying Td,i model are quite different than the constant 
T o j model; 2) the quality of the fits improve to the point 
of be comparable to halo fits; and 3) the predicted T^j 
variations agree qualitatively with observed variations. 
It is possible that a more quantitative match could be 
recovered with the proper observations. 

In conclusion, while we find little statistical evidence to 
support the claim that Milgrom models explain the be- 
havior of this sample of galaxies more economically than 
do dark matter models, there are several sources of sys- 
tematic uncertainty that complicate the situation. Given 
these systematic uncertainties (distance, gas content, T 
gradients), it is difficult to see how Milgrom mass mod- 
els are of lower quality than dark matter halo models. 
However, if one believes that stellar population models 
accurately describe galaxies, then correct models of the 
stellar components of real galaxies should reveal trends 
expected from such models. This is one piece of evidence 
in favor of Milgrom over dark matter models; Milgrom 
model T values show more of a trend with color than do 
dark matter model T values for the same galaxies. 

Dynamical analyses which purport to rule out the Mil- 
grom relation must be done on carefully selected galaxies 
for which systematic errors in dynamical assumptions are 
minimized. It is easy to demonstrate superior dynamical 
fits by dark matter halos than by the Milgrom relation, 
but on closer examination it is clear that in the majority 
of cases, the discrepancy can be explained by reasonable 
assumptions about particular known systematic errors. 
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APPENDIX 
BAYESIAN STATISTICS PRIMER 

The application of Bayesian reasoning to statistical problems in astronomy has become much more widespread 
in recent years, and numerous textbook treatments exist. Here we give only a brief summary of Bayesian model 
comparison used in this paper, and point the interested reader to I Javn es ( 2003) ; iBretthorstl (|1990l ). 

Bayesian statistics is based on Bayes' theorem, 

P(M\DI) is read as "the probability that M is true given D and I". Let us make this concrete. Say that proposition 
M is a given model, proposition D is the data fit by the model, and proposition / is any initial information we have 
about the system. The left hand side of Equation IA1I is then the probability that a model is true given the data and 
prior information. If we have one data set, two models, and the same prior information, we can write Eg uat ion I A 1 1 for 
each model Mj and M k . Since D and / are the same in both equations, we can divide the two to cancel P(D\I). We 
now have a ratio of the probabilities of the models being correct, 

PjMjlDI) _ PjM^PjDWjT) 

P(M k \DI) P(M k \I)P(D\M k iy [ > 

Following I Javnesl (|2003l ) we have that, 

P(D\MI)= J d®P(D\@MI)P(@\MI) 

= J d®L(®)P(@\MI), (A3) 

where are the (1...N) free parameters, L is the likelihood of a model, and P(®\MI) is the prior probability. If we 
normalize the likelihood by its maximum value L max , we have that P(D\MI) = L max W , where W is the Ockham 
factor, 

W = f d®^^-P(@\MI). (A4) 



In order to calculate the Ockham factor we need to decide what to use for the prior parameter probability. In the 
absence of any specific knowledge of the parameters, we choose a top-hat prior, i.e., a constant value p in a specified 
range. Since the prior probability must be normalized, we can easily calculate what p is. J P(®\M I)d® = 1 implies 
that p = 1/(A0iA02...A0tv), where A0j is the range of a given free parameter. With this prior we can write, 

P(D\MI) = L max p f d®^l 

= L maxP V(Q), (A5) 
where V(£l) is just the integral of the normalized likelihood over all parameter space and is related to the volume of 
allowed parameter space. We are now ready to rewrite the right-hand side of Equation IA2( 

P(M 3 \DI) = P(Mj\I) L maXij Pj VjjQj) 

P(M k \DI) P{M k \I)L m ^ kPk V k (n k )- [ ' 

Assuming our prior information leaves us equally ignorant about both models, the first term on the right-hand side 
of Equation I A6I is 1. We are left with an expression that is simple and sensible. If Mj provides a better fit (i.e., has 
a smaller xV) than M k , then the likelihood ratio (the second term) is large, thus pushing the ratio in favor of Mj, 
as one would expect. If the two models have the same free parameters, then the third ratio relating the parameter 
ranges is 1 and the last ratio is w 1, We are left with the obvious answer that the model with the highest likelihood is 
the more probable model. However, the real power of this approach emerges when you compare models with different 
numbers of parameters. For concreteness, let us say that Mj has two free parameters a and (3 and model M k has two 
additional parameters 7 and 8. We know L max for each model, and we assume that L max .fc < L max ,j- We also know 
that, pj = l/(AaA/3) and p k = l/(AaAf3A^A5), so the p ratio term is A'jAS. Since M k has the larger number of 
parameters, let us suppose that Vj(Qj) > V k (ft k ). The likelihood ratio is < 1, the V ratio is > 1, and the p ratio is 
dependent on the ranges that we believe two of the parameters fall in. The Bayesian model probability ratio does not 
necessarily choose the model with the highest likelihood: it determines if the benefit of adding more fitting parameters 
balances the cost of making the model more complicated. 
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TABLE 1 

Best-fit \\ values, T values (with statistical uncertainties a), and systematic uncertainties (s) for Milgrom models. The 

LAST COLUMN CONTAINS THE RATIO BETWEEN THE KINEMATIC DATA'S RADIAL EXTENT i?k,max AND THE DISK SCALELENGTH Td- 



Galaxy 



^6 



Tj),. 



ESO 215g39 


1.48 


1 


ESO 216g20 


4.16 


2 


ESO 263gl4 


1.93 





ESO 267g29 


5.42 


2 


ESO 268g37 


4.11 


1 


ESO 268g44 


2.05 


1 


ESO 317g41 


4.58 


2 


ESO 322g36 


4.53 





ESO 322g42 


3.67 





ESO 322g44 


1.61 





ESO 322g45 


2.16 


1 


ESO 322g76 


2.05 


1 


ESO 322g77 


3.25 


3 


ESO 322g82 


4.53 


1 


ESO 323g25 


2.81 


2 


ESO 323g27 


2.83 


1 


ESO 323g39 


4.55 





ESO 323g42 


1.66 





ESO 323g73 


3.10 


1 


ESO 374g03 


3.29 





ESO 375g02 


2.08 


1 


ESO 381g05 


2.75 


5 


ESO 382g06 


2.48 


3 


ESO 435g26 


11.98 


1 


ESO 437g04 


2.91 


2 


ESO 437g31 


1.07 


2 


ESO 438g08 


2.32 





ESO 438gl5 


2.60 





ESO 439gl8 


4.14 


2 


ESO 439g20 


1.86 


2 


ESO 444g47 


1.69 


2 


ESO 445gl5 


3.03 


2 


ESO 445gl9 


2.31 


1 


ESO 445g39 


4.22 


5 


ESO 446g01 


2.23 


1 


ESO 501g01 


2.33 


1 


ESO 501g68 


2.13 


2 


ESO 502g02 


1.74 


2 


ESO 509g80 


3.75 


2 


ESO 569gl7 


1.67 


2 



8:81 



14 



8:8? 
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8:?§ 
8:&g 



8:8i 



Mi 



Mi 



8:81 

8:?e 

8:41 



Mi 

i 
S 

8:8? 
8:83 



0.63 
0.14 
0.08 
1.77 
0.64 
0.08 
0.04 
0.47 
0.02 
0.06 
0.02 
0.27 
0.13 
0.17 
0.83 
0.13 
0.84 
0.01 
0.29 
0.04 
0.24 
2.86 
1.11 
1.05 
0.23 
1.49 
0.29 
0.10 
0.37 
0.17 
0.15 
0.17 
0.03 
0.63 
0.50 
0.44 
0.31 
0.02 
0.19 
1.83 



2.79 
1 31 



2 14 04 
z - lqi 0.02 



3.09 
1.43 
1.09 



0.27 

8:?! 
8:A8 



0K 



10 



0.72 1 
1 43 

m 

2.08 
7.99 g 
15.56 ' 



i 4 i, 

1.26 

5 10 26 
°- iU 0.31 



1 60 28 
l.OU a21 



1 64 32 
L -°^ 0.36 



3.42 
3.20 
2.42 
2.07 
3.66 
1.66 
2.20 
2.90 
2.00 
1.58 
1.12 
2.89 

1.32 
4.76 
2.08 
3.17 



1.61 



■■■m 



m 

i.36 
0.03 

Mi 
8:81 
8:9§ 



S ?B.I 

0.57 
0.37 

1.90 

0.15 
0.04 
0.38 

2.33 
0.02 
0.19 
0.23 
1.72 
8.63 
0.05 

0.08 



0.43 

0.72 
1.38 
0.23 
1.65 
1.53 
0.09 
0.37 
0.42 
0.15 
0.13 
0.08 
0.76 

0.44 
0.15 
0.02 
0.14 



.max 

3.35 
5.20 
4.64 
3.81 
3.78 
4.45 
5.09 
3.71 
3.66 
4.05 
2.11 
4.47 
6.33 
4.12 
4.75 
4.25 
3.09 
4.50 
5.11 
3.84 
5.35 
3.73 
4.03 
6.04 
4.51 
3.74 
4.71 
6.62 
4.46 
3.98 
4.95 
5.17 
3.69 
4.11 
3.97 
2.76 
3.84 
4.50 
4.83 
5.72 



a T values have units (Mq/Lq). 
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TABLE 2 

Logarithms of probability ratios for the various models. The columns list the relative probabilities of the the various 
models to mllgrom models. negative values indicate that the mllgrom model is preferred. the largest positive value for 
a galaxy comes from the model (other than mllgrom) that best describes the data. 



Galaxy 


-^stars — only 


A PI 


A NFW 


Apowcr — law 


ESO 215g39 


-44.66 


40.26 


32.27 


23.78 


ESO 216g20 


-32.43 


9.06 


3.92 


5.61 


ESO 263gl4 


-49.46 


7.16 


4.38 


0.29 


ESO 267g29 


-47.16 


27.47 


24.76 


25.42 


ESO 268g37 


-170.16 


33.57 


17.81 


1.97 


ESO 268g44 


3.21 


8.95 


6.34 


0.45 


EbO 317g41 


0.05 


4.24 


-15.56 


-17.21 


ESO 322g36 


-144.46 


16.31 


15.70 


14.20 


T ' { ' /" \ Oil,'! 1 i"> 

EbO 322g42 


-81.50 


49.17 


11.56 


12.90 


ESO 322g44 


-153.87 


7.76 


3.93 


5.06 


EbO 322g45 


-39.33 


14.59 


11.65 


13.18 


ESO 322g7b 


-17.42 


0.13 


-4.39 


-1.50 


ESO 322g77 


-3.96 


8.18 


7.07 


8.02 


ESO 322g82 


158.84 


162.29 


163.86 


158.06 


ESO 323g25 


59.15 


177.19 


195.25 


113.44 


ESO 323g27 


-4.76 


21.32 


17.77 


12.86 


ESO 323g39 


-55.35 


20.97 


12.84 


11.75 


ESO 323g42 


36.22 


79.43 


64.07 


61.33 


ESO 323g73 


-94.97 


53.64 


51.04 


52.33 


ESO 374g03 


-78.04 


71.89 


30.20 


0.59 


ESO 375g02 


-84.92 


21.97 


21.95 


20.80 


ESO 381g05 


-128.24 


215.20 


216.79 


212.01 


ESO 382g06 


10.70 


9.89 


6.46 


7.60 


ESO 435g26 


126.10 


182.98 


192.48 


126.17 


ESO 437g04 


-153.84 


207.50 


171.75 


158.56 


ESO 437g31 


-67.01 


13.52 


-3.45 


-2.09 


ESO 438gU8 


77.78 


76.00 


73.35 


74.73 


ESO 438gl5 


-92.73 


43.05 


39.63 


41.19 


ESO 439gl8 


-114.20 


73.18 


63.88 


35.15 


ESO 439g20 


-65.38 


35.96 


24.74 


3.19 


ESO 444g47 


-96.68 


14.98 


2.95 


5.09 


ESO 445gl5 


-79.48 


-9.18 


-13.31 


-13.43 


ESO 445gl9 


-67.36 


7.28 


13.19 


13.93 


ESO 445g39 


-0.26 


7.52 


0.03 


1.37 


ESO 446g01 


-219.48 


4.36 


0.25 


1.24 


ESO 501g01 


-66.08 


31.50 


8.88 


-11.19 


ESO 501g68 


-4.56 


-1.74 


-4.70 


-2.35 


ESO 502g02 


-80.73 


16.47 


12.31 


5.60 


ESO 509g80 


56.56 


169.81 


166.45 


90.89 


ESO 569gl7 


5.67 


43.13 


46.70 


9.39 
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R (kpc) R (kpc) 

Fig. 1. — Comparisons between Milgrom relation and pseudo-isothermal halo rotation curves. The solid lines are the Milgrom rotation 
curves and the dot-dashed li nes are the pseudo-is othermal halo curves. The error bars mark the circular speeds derived from the velocity 
map during kinematic fitting (Barnes & Sellwood 2003). Panels a) (ESO 322g76) and b) (ESO 445gl5) illustrate cases in which the Milgrom 
fit is preferred over or comparable to the other models. Panels c) (ESO 268g44) and d) (ESO 438g08) show two cases where the Milgrom 
models are not favored over any of the other models. ESO 268g44 has substantial spiral structure and ESO 438g08 appears to have an 
interacting companion. Panels e) (ESO 323g42) and f) (ESO 382g06) are the axisymmetric galaxies for which the Milgrom fit is never 
favored. 
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Fig. 2. — The color-To / correlations for the galaxies in this sample assuming Newt onian star s -only models (a) and Milgrom models 
(b). The dashed line is the best linear fit to the points and the solid line is the fit from [Bell et al\ J2003I) . The error bar is the estimated 
uncertainty in galaxy color. 
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Fig. 3. — Th e color-Tp , / corr elations for the Hydra cluster galaxies. Again, the dashed line is the best linear fit to the points and the 
solid line is the lBell et "oil '(20031 fit. 
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Fig. 4. — Comparisons between Milgrom (solid), pseudo-isothermal (dot-dashed), and Newtonian stars-only (dotted) rotation curves for 
ESO 3 23g42 (a) and ESQ 382g06 (b). The points with error bars show the circular speeds derived from the velocity map during kinematic 
fitting (Barnes & Sellwood 2003). 
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Fig. 5. — Comparisons between spatially-constant Trj r Milgrom (solid lines), pseudo-isothermal (dot-dashed lines), and variable Trj j 
Milgrom (triple-dot-dashed lines) rotation curves for ESO 322g76 (a), 323g42 (b), and ESO 382g06 (c). Again, the points with error bars 
denote the circular speeds derived from kinematic fitting of the velocity map. 



